Isospin excitations of a trapped ID gas of attractively interacting fermions 
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We consider a gas of fermions with a short-range attractive intercomponent interaction in a 
parabolic external potential and derive the conditions of the local density approximation. The 
obtained spectrum of quasiparticle (isospin) excitations shows equidistant low-energy levels, which 
is equivalent to a linear momentum dependence and is fundamentally different from the ordinary 
Dirac spectrum in the spatially uniform case. 
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Fast progress in experiments with cold atoms has led 
to the creation of one-dimensional (ID) atomic gases by 
(tightly) confining the motion of atoms in two direc- 
tions to zero point oscillations. One-dimensional quan- 
tum gases show a remarkable physics not encountered 
in higher dimensions. In particular, since the density 
of states in ID increases towards the zone boundary, 
the effective strength of interactions in the ID gas in- 
creases with decreasing density. Therefore, decreasing 
density in a gas of attractively interacting fermions one 
can crossover from the BCS-like regime of strongly over- 
lapped pairs to the regime of compact pairs forming a 
Bose-Einstein condensate (EEC). For spatially uniform 
systems this problem is well understood due to avail- 
ability of exact solutions obtained by the Bethe Ansatz 
in combination with powerful bosonization tech- 
niques [J, 5, 6]. 

The ID gases are usually obtained in an external har- 
monic potential, which introduces a finite size of the sys- 
tem and makes it spatially nonuniform. Trapped ID 
Bose gases have been intensively studied in the last years 
(see [7| for review). Recently, the ID regime has been 
achieved for atomic fermions and the discussion of 
trapped ID Fermi gases was focused on the occurrence 
of the BCS-BEC crossover, manifestation of spin-charge 
separation, role of the imbalance between atomic compo- 
nents (see for review). One of the key problems is 
revealing the influence of an external harmonic potential 
on the quasiparticle spectrum of this system. 

In this paper we consider a ID TV-component Fermi 
gas with a point-like attraction in a parabolic trap. We 
first treat this problem in the limit of iV >> 1 using the 
standard 1 /TV-expansion ll|. In the limit of — > oo 



the saddle point approximation becomes exact and one 
can deal with the spatially non-uniform distribution of 
particles in a controllable way. We show that the spatial 
inhomogeneity strongly affects the quantization rules for 
the quasiparticle spectrum. As a result, the system in 
a parabolic trap cannot be mapped onto an integrable 
system in a rectangular box. We then discuss the ap- 
plicability of the large- A^ results for the experimentally 



relevant case of a two-component Fermi gas (A^ = 2). 

In terms of field operators ^jJ^ and ^pj , with the index 
labelling the fermionic species, the bare Hamiltonian is 



H = ^ j dx 



2m 



2 9 

muj X 



(1) 



3' 1^3 



where m is the atom mass, lo is the trap frequency, the 
coupling constant of the point-like attraction is written 
as —g/N with g > 0, and we set h ~ 1. We will work in 
the thermodynamic limit defined by the relations 



oAf /N = const, TV — > cx), w — > 0; 



(2) 



with A/" being the total number of particles, and Ep the 
Fermi energy. 

We first briefly outline the results for the uniform case 
where the system is in a rectangular box with periodic 
boundary conditions. Then, in the thermodynamic limit 
the excitation spectrum consists of a gapless collective 
mode describing fluctuations of the total density (charge 
mode) and A^ — 1 branches of particles (isospin modes) 
[l^. [isj. The latters have spectral gaps: the particle of 
the q-th branch transforms according to the fundamen- 
tal representation of the SU(N) group with the Young 
tableau consisting of one column of height q. In the limit 
of 7 << I where 
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and n is the mean density of each fermionic component, 
the spectrum of the gaped particles is approximately 
Lorentz invariant: 



e,{p) ^ ^{pvpf+Ml (4) 

Afg = Mosin(27r(j/Ar), q^l,...N-l; 

Mo = CEp-f'^'' exp(-2^/7), (5) 

where C is a numerical constant. At a finite A^ these 
particles interact and the ones with q > 1 are bound 
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states of the fundamental particle with q = 1- However, 
in the limit of ^ oo the interaction vanishes and the 
only particles which remain are the particle with q = 1 
and its antiparticle having q = N —1. The case of iV >> 



1 can be treated perturbatively using 1 /N expansion 11 1 . 

This approach can be employed for a spatially nonuni- 
form system as is done in this paper. To make calcula- 
tions easier we assume that in the major part of the trap 
the coupling constant ([3|) in which n now depends on the 
coordinate, is small. This allows us to use the Thomas- 
Fermi density profile for non interacting fermions: 
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where xq is the Thomas-Fermi (half)size of the sam- 
ple. The Thomas-Fermi approximation breaks down 
at distances near the Thomas-Fermi boundary, where 
dn-\x)/dx - 1 and, hence, (1 - x^/x^)^^^ - Af-^^^. 
The condition of weak coupling breaks down in a narrow 
range of distances near the boundary, where "/{n{x)) ~ 1 
and (1 - x^/a;^)^/^ ~ 7(0) = mg/nno < 1. 

Assuming that the weak coupling approximation holds, 
we can linearize the spectrum introducing right and left 
moving fermions R, L: 

^(;j{x) = e^^'=^(^)^i?j(.T) +ci'=^(^)^i:j-(x). 
The coordinate-dependent Fermi momentum is kp{x) = 



Tm{x), and Hamiltonian ((T|) takes the form: 

H = ^ r drUi-R+drRj + L+drLj] 
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where 7 is given by Eq. and 

dx' /vp{x') = u!^^ arcsin(a;/a;o), 



(8) 



with the Fermi velocity vp — kp/m, and tq = 7r/2a;. 

At >> 1 the interaction can be treated in the 
saddle-point approximation which becomes exact in the 
limit of iV — > cxD,7 = const. To develop such an ap- 
proximation one does the Hubbard-Stratonovich trans- 
formation introducing an auxiliary complex field A(t, t), 
and formally integrates over the fermions. The resulting 
effective action is 



I 27(t) 
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i{dt ^dr ) A 
A* m + dr) 



where t is the Matsubara time. Expanding the logarithm 
in gradients of A we obtain the Lagrangian density: 
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where A ~ Ep is the high-energy cut-off, and the dots 
stand for the terms containing higher powers of da^^^ 
(a = t,T). These terms are small under the condi- 
tion |Ao| = |A(0)| ^ cj, since the characteristic scale 
of T is . This is equivalent to having the system 
size xq much larger than the largest correlation length 
Co = «f(0)/|Ao|. 

We now write A in the form 



pexp(i(? 
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and obtain the action in which fluctuations of the phase 
(j) are decoupled from the fluctuations of the amplitude 
|A|^p: 



c = m - 

N 
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7(t) 



- ln(A/p) 
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P i-' 

There are two features of this expansion, with a different 
level of robustness under deviations from the condition 
of weak coupling, 7 << 1. First of all, the phase (j) 
is decoupled from the amplitude field, which is a sim- 
ple consequence of (p being the Goldstone mode of the 
Charge Density Wave order parameter field ([TO)) . As such 
it should be gapless and weakly coupled to other excita- 
tions. However, the second feature, namely the fact that 
the stiffness of the field ip is independent of the parti- 
cle density n{x), holds only in the limit of 7 << 1. As 
follows from the Bethe Ansatz calculations valid in the 
uniform case, the stiffness (or the Luttinger parameter 
Kc) starts to acquire the density dependence at 7 ^ 1. 

Since the action following from Eqs. pT jl -lfTS ]) is pro- 
portional to large A'', fluctuations are suppressed. In par- 
ticular, one can find excitation energies by considering 
the Dirac Hamiltonian with the coordinate-dependent 
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mass determined by the minimum of the action (|13p : 

2p[27r/7(T)-hi(A/p)]-p-i9^^hip+... = (14) 
The sohition can be represented as 

p{t) ^ p{0) exp \ 2TT'^[n{0)~n{T)]/ gm+^drr \np + . 



2p2 

The second term in the exponent of this expression can 
be omitted if one satisfies the inequaHty 

Ao/uj > 1/^, (15) 

where 70 = 7(0). One then has 

p{t) = Aoexp{27r2[no -n{T)]/mg}. (16) 

Now we find the excitation spectrum of quasiparticles. 
As we have said, in the hmit of — > 00 their spectrum is 
decoupled from the phase excitations. Then the fermions 
are described by the effective Dirac Hamihonian 



Rj 



,(17) 



where p(t) is given by Eq. (|16p . A uniform system in a 
rectangular box with periodic boundary conditions, is de- 
scribed by the same Hamiltonian, with p(t) — Aq. Equa- 
tions of motion following from this Hamiltonian yield so- 
lutions in the form of plane waves with the spectrum 
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[■Kkvp/lf + Al, fc = 0,±l,±2. 



(18) 



where 21 is the length of the box. This corresponds to the 
spectrum of g = 1, A^ — 1 excitation branches in Eq.Q. 

Actually, having derived Eq. for pir) entering 

the Hamiltonian Hp ([TT]), we obtained the local density 
approximation for our nonuniform system. Our results 
show that in the Dirac Hamiltonian p?)) one can use the 
same expression for the mass (gap) as in the uniform 
case, but the exponent is coordinate dependent through 
the spatial dependence of the density (parameter 7) and 
the coordinate dependence of the preexponential factor 
is omitted. However, the quasiparticle spectrum is quite 
different from that of Eq. p^ . 

Eigenfunctions of Hamiltonian (jl7p can be written as 



R = p^^^F, L={eR + iRr)lp, 



(19) 



p'{r)-,^+(P^-(P^] +i.^) F = 0,(20) 



where the function F{t) satisfies the equation 

2 

r „ „ / n / n_ \ 

' ' ■ ' ' 2p \2p 

and the notations F^. and Ft^ mean the first and sec- 
ond derivative with respect to r. For finding the lowest 
eigenstates one can approximate 



p\r) « 



1 
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and check that under the condition p5|) the terms in 
the round brackets in Eq. (|20|) are small at least as 
(tii/^yTlJAo)^/^ compared to both terms on the rhs of 
Eq. (PT|) . Then, omitting the terms in the round brackets, 
for the energy eigenvalues we obtain: 



tu = ^A2 + Aoc^(8V7o)i/2(fc + 1/2), (22) 

where fc > is an integer. Identifying fc as a (rescaled) 
momentum of a particle confined in a rectangular box 
and comparing Eq. with Eq. one sees that the 
quasiparticle energy levels in the parabolic trap cannot 
be mapped onto the Dirac particle spectrum in the box. 

For finding high-energy eigenstates one can use the 
WKB quantization rule. Omitting again the terms in 
the round brackets in Eq. (120]) we have 



dT\/ — p^ij) = Trk, 



(23) 



where T(e) is determined by the condition e — p{t). 
Rewriting Eq. ITBl) as 

p(t) — Aq exp{27r(l — cosijjr)/7o} 

we notice that considering e <ti Ep one can write the 
exponent of p(r) as -kllP-t^ j^t^. This leads to 



r(e)«-J-ln(e/Ao). 

V TT 



Then performing the integration in Eq. 
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we obtain: 



= 7rfc, (25) 



where for eigenstates near the bottom of the gap, that is 
for (efc — Aq) ^ Aq, the coefficient /3 is equal to 7r^/16 
and Eq. ((25)) gives the same result as Eq. ((22l) at large 



quantum numbers k: ef. 



2 - Ag -f Aot^(87r/7o)i/2fc. For 
excitation energies ek ^ Aq one has f3 close to unity, 
and in the limit of 3> Aq corresponding to quantum 
numbers k 3> Aq^/jo/^^, Eq. (^5)1 yields: 
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470/ VMR727oF^^^fc7Ao] 



(26) 



As we see, this result also differs significantly from the 
simple Dirac spectrum. 

In Fig. 1 we compare the results of Eqs. (|^ and P5|) 
with numerical calculation of the spectrum from Eq. (|20p . 
One sees a good agreement even for Aq/uj about 20, and 
for larger values of Aq/w the agreement is nearly perfect. 

We believe that our main result, a pronounced dif- 
ference between quantization rules for the quasiparticle 
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FIG. 1: Energy levels for: (a) Aq/cj = 26 (70 = 0.44); (b) 
Ao/o; = 95 (70 = 0.5). The solid and dashed curves show the 
results of Eq. ([22]) and the results of Eq. ^ with f3 = 1. The 
crosses are the results of numerical solution of Eq. (|20|) . 



spin-up and spin-down particles get kicks in the oppo- 
site directions, which creates isospin modes. At the same 
time, charge (density) modes corresponding to in-phase 
oscillations of the two components are not excited. 

In conclusion, we have found the quasiparticle (isospin) 
spectrum for attractively interacting fermions in a 
parabolic potential. The spectrum shows equidistant 
low-energy levels (linear momentum dependence) and is 
drastically different from the ordinary Dirac spectrum in 
the spatially uniform case. Experimental verification of 
this result will provide a clear demonstration of the fact 
that the parabolic confinement can fundamentally change 
the properties of the quantum gas. 
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spectrum in harmonic and rectangular confining poten- 
tials, remains valid in the experimentally relevant case of 
a two-component Fermi gas {N — 2). It is likely that the 
local density approximation for the mass gap p calculated 
along the lines of 14, 1^, together with Eq. ((23)) . remain 
robust. The main difference of small from large TV is a 
dramatic increase of interactions between quasiparticles. 
This, however, will not have a strong effect on the energy 
levels as soon as the number of excited quasiparticles is 
macroscopically small. The interactions affect commuta- 
tion relations of the quasiparticle creation and annihila- 
tion operators (see, e.g.. Chapter 34 of 5] and references 
therein) and this leads to changes in the behaviour of 
matrix elements and correlation functions [l^ . 

For a trapped ID atomic ultracold Fermi gas one can 
expect the number of particles J\f ^ 10"^. Then, re- 
calling that in the two-species gas the Fermi energy is 
Ef = Afuj/2, for 70 « 1 we obtain Aq/lu « 10 and our 
results are applicable. The isospin modes can be excited 
optically, for example by using pulses of polarised 
light acting on one of the atomic components (spin-up) 
and pulses of light acting on the other component 
(spin-down). The a~ and cr^ pulses provide periodic op- 
tical potentials shifted by a quarter of a wavelength with 
respect to each other, so that the minimum of the po- 
tential corresponds to the maximum of the cr"*" potential 
17| and the sum of the two potentials is zero. Thus, the 
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